Edge solitons of topological insulators and fractionalized quasiparticles in two 

dimensions 
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An important characteristic of topological band insulators is the necessary presence of in-gap edge 
states on the sample boundary. We utilize this fact to show that when the boundary is reconnected 
with a twist, there is always zero-energy defect states. This provides a natural connection between 
novel defects in the two-dimensional p x + ip y superconductor, the Kitaev model, the fractional 
quantum Hall effect, and the one-dimensional domain wall of polyacetylene. 
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Excitations carrying fractional quantum numbers (e.g. 
fractional charge), such as the quasiparticles in the frac- 
tional quantum Hall effect [l[ , have always been an sub- 
ject of interest. In 1976 Jackiw and Rebbi Q wrote a 
seminal paper which laid the foundation of charge frac- 
tionalization in one spatial dimension. Four years later 
the influential paper of Su, Schrieffer and Heeger Q pro- 
poses the "Jackiw-Rebbi soliton" as the charge carrier in 
doped polyacetylene. Today, all quantum number frac- 
tionalization phenomena, such as the fractionalization of 
magnon into spinons Q , can be attributed to the Jackiw- 
Rcbbi mechanism. 

In the last fifteen years, starting with the fractional 
quantum Hall effect, condensed matter physicists stum- 
bled upon several instances where quantum number frac- 
tionalization occurs in two spatial dimensions. These in- 
clude the quasiparticles of the "Pfaffian" quantum Hall 
state 0], the vortices of a spin-polarized p x + ip y su- 
perconductor [f| 0, H, 0, EH > and the topological excita- 
tions in a spin model proposed by Kitaev [ll|. However, 
what is lacking is a general framework, like the Jackiw- 
Rebbi theory in one dimension, specifying the condition 
under which fractionalized excitations will appear. In 
this paper we provide such a mechanism and reveal its 
connection to the Jackiw-Rebbi theory. In particular, we 
show that, in two dimensions, fractional charge will nat- 
urally appear around defects of "topological" insulators 

[3 El m- 

Recently it has been shown that the existence of frac- 
tional charge in the quantum Hall effect is connected to 
the existence of fractionally charged domain wall in cer- 
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In the 



tain one dimensional systems 
following we generalize such connection and show that 
when the boundary of a topological insulator, an insu- 
lator which necessarily possesses in-gap edge state, are 
reconnected with a twist, there is always zero-energy de- 
fect states possessing fractional quantum number. 
Let us begin by considering the Kitaev model 



This exactly soluble model describes a honeycomb lat- 
tice of quantum one-half spins interacting via three type 
of nearest-neighbor interactions (Fig.QJ,). The Hamilto- 
nian is given by 
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where "w/b" abbreviates for the white/black sublattice 
(see Fig.QJ,), and n+e M is the nearest neighbors of n along 
the fj, bond, are the Pauli matrices. By perfor ming 
Jordan- Wigner transformation it was shown in Ref. [191 ] 
that this model is equivalent to a free Majorana fermion 
model: 

+ ^3-Dn7n+e 3 7n}; ( 2 ) 

new ^l— 1,2 

where D n — ±1 is a classical Ising variable and 7 n 's are 
Majorana fermion operators [19J. Since the honeycomb 
lattice is consisted of two sublattices, Eq. can be re- 
cast into 
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where — (7,1 , 7n+e3 

). In Eq. J3) the 2 x 2 coupling 
matrix H nm has value: H nm — J^D n <j2 for n = m, and 
Hum, = z f('iJfi/2)< J ± for n ^ m. The upper sign applies 
if the bond linking m to n is a black-to-white bond, and 
the lower sign applies if it is a white-to-black bond. 

In the ground state, the classical Ising variable takes 
value D n — 1 modulo a global flip per row 19]. For D n = 
1, Hm is translation invariant, and can be diagonalized 
by Fourier transformation. The Bloch matrix of Hm is 
given by HmQ*.) = hi(k)ai + /i2(k)(T2 with 

ft-i(k) = — J 2 sina(k) + Ji sin/?( k), 

h 2 (k) = J 3 + J 2 cosa(k) + JiCos/3(k), (4) 



FIG. 1: (a) A graphical representation of the Kitaev model. 
There are two sublattices (white and black), and three types 
of bonds (labeled by 1,2,3, or black, red, blue), (b) The graph- 
ical representation of H+Ht- In the Majorana representation, 
the three-spin interactions in Eq. ([5]) become next neighbor 
hopping along the zigzag chain in the x-direction. The arrow 
represents the direction in which the second neighbor hopping 
matrix elements are —1J4. 



and a(k) = (y/3k x - 3k y )/2, (3(k) = (V3k x + 3k y )/2. 
It is easy to show that for J\ + J 2 > J3 > \J\ — Ji\ 
the energy spectrum is gapless. The possible connection 
to the quantum spin liquids has been discussed in Ref. 
[20I ] . In Ref. 11 1 a magnetic field is introduced to open 



an excitation gap in this parameter region, and the non- 
Abelian quasiparticles become low energy excitations of 
this gapped phase. 

Unfortunately the magnetic field spoils the integrabil- 
ity of the model. Here we propose a different way of 
opening a gap while maintaining the integrability. This 
is achieved by adding the following three-spin interaction 
to Eq. Q 
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FIG. 2: The /i-surface for the Kitaev model (a) and the 
Px + ip y superconductor (e). In panels (b) and (f) the h- 
surfaces are dissected to expose the origin (red dot) . (c) Topo- 
logical defects in the Kitaev model. To the right of the grey- 
shaded plaquette the sign of D n 's are reversed. As a result the 
corresponding vertical hopping matrix elements change sign. 
They are shown by the dashed bonds, (d) The eigenspectrum 
assoicated with two far separated topological defects in the 
Kitaev model with Ji = J2 = 1, J3 = 0.2 and J4 = 0.5. (g) 
A vortex centered at the gray-shaded plaquette in a p x + ip y 
superconductor. The pairing order parameter is shown for 
four bonds, (h) The energy spectrum of two far-separated 
vortex-antivortex pair with t = 1, Ao = 0.5 and fi = 0.3. 



For a /i-surface not containing the origin, there is an 
integer topological index 



1 

8^ 
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Here (ijk) denote three adjacent sites (with i being the 
left-most one) along the zigzag chain running along the 
x-direction. Depending on whether (ijk) form an up- 
pointing or a down-pointing triangle, we use the first or 
second term of Eq. ([5]). In terms of the Majorana fermion 
operators this amounts to adding a second nearest neigh- 
bor hopping between sites along the zigzag chain: 
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Then the Bloch matrix becomes ifjf(k) = /ii(k)<7i + 
h 2 (k)a2 + h 3 (k)a 3 , where h 3 (k) = 2 J4 sin(v / 3fc K )- The 
vector function h(k) is a continuous mapping from 
the first Brillouin zone to the space spanned by h = 
{hi, h2,h 3 ). The image is a closed two-dimensional man- 
ifold (henceforth referred as the ft,-surface). Since the 
eigenvalues of the Block matrix are ±|h(k)|, it follows 
that if the /i-surface contains the origin, the spectrum is 
gapless, otherwise the spectrum has a gap. 



which counts the number of times the unit vector h wraps 
around the origin. As shown in Ref. [2l|. V is propor- 
tional to the well known "TKNN" index [22| in the case 
of two bands. Spectra characterized by different V are 
topologically distinct. They can not be deformed into 
each other without gap closing. In the parameter regime 
where non-abelian quasiparticle exists V = 1. 

Topological excitations of the Kitaev model are cre- 
ated by reversing the sign of D n 's in Hm along half a 
row. This is shown by the dashed bonds in Fig. J2fc). In 
Fig. ((2ji) we have shown the result of numerical diago- 
nalization for a system of 1600 sites with toric boundary 
condition. Because of the boundary condition two de- 
fects are introduced, they are separated by 100 sites in 
the x-direction. They introduce two mid-gap states with 
a tunnel-splitting (which is already invisible here) which 
decreases exponentially with the separation. 

Now we switch gear to discuss the vortices in a 
spin-polarized p x + ip y superconductor. Let us consider 
this problem on a square lattice. Similar to the Kitaev 



model, the Hamiltonian can also be written in the form 



and 



is a fermion 



of Eq. Q, except = (c+, 
operator. The Bloch matrix is characterized by h(k) = 
(— A smfc y , A sin k x , —t(cosk x +cosk y )~ii). Here 
A is the pairing amplitude, t is the hopping integral, 
and fi is the chemical potential. The /i-surface is shown 
in Fig. ([2b-) and (J2f ) . Straightforward calculation shows 
that V = 1. After a singular gauge transformation, a 
vortex can be created by reversing the sign of the hop- 
ping matrix elements along a cut as shown in Fig. ([2fe) . 
Explicit calculation shows that there is also a zero mode 
associated with each vortex (Fig.[2fr)- 



For each free Majorana or Bogoliubov fermion model 
discussed above, there is a free fermion model with iden- 
tical excitation spectrum. To obtain this fermion model 
we simply replace 4^ in Eq. ([3|) by the fermion operator 
= ( c in' c 2n) where 1 and 2 are "flavor" indices (they 
might represent the two sites in the unit cell of a lattice). 
This fermion model acts as a representative of all models 
which share the same H nm (hence the same eigen spec- 
trum). However, while the representative fermion model 
is global U(l) invariant, the Majorana and Bogoliubov 
fermion models only have Zi symmetry. The fact that 
in the Z2 models the particle number is only conserved 
modulo two is the root of non-abelian statistics. In the 
rest of the discussion, we refer to the system described 
by a gapped free fermion model with non-zero P as a 
"topological band insulator". Thus the representative 
fermion models for the Kitaev model and the p x + ip y 
superconductor are topological band insulators. Know- 
ing the properties of edge states and defects of the rep- 
resentative fermion model, one can readily deduce the 
corresponding properties of the Majorana fermion (Ki- 
taev) or the Bogoliubov fermion (p x + ip y ) models with 
the same H nm . For example, while in the fermion model 
the edge states are free fermions, and the defect zero 
modes carry half fermion quantum numbers, those in the 
Majorana/Bogoliubov fermion model are true Majorana 
fermions. 



In the following we provide an unifying mechanism for 
the appearance of defect zero mode when the representa- 
tive fermion model describes a topological band insulator. 
As an example, let us consider the fermion representative 
of the Kitaev model (Fig.[3K). Fig. ([3b) shows the gapped 
spectrum as a function of momentum along the longitu- 
dinal circle. Due to its topological nature, if we remove 
a row of bonds (Fig. [3b), in-gap edge states appear [23] 
as shown in Fig. ([3b!). The left and right moving chiral 
edge fermions, represented by the arrowed red lines in 
Fig. ([3b) , are described by a massless Dirac Hamiltonian 
in ID. If we reconnect the two edges, but with weaker 
bonds, a smaller gap reappear in the edge spectrum (see 
Fig. [3b and Fig.[3f). The edge fermions are now described 





FIG. 3: (a) The Kitaev model and (b) its eigen-spectrum 
as a function of wavevector k x along the x-direction around 
k x — tv. (c) One row of bonds along x are removed and (d) the 
corresponding energy spectrum. The arrowed-red lines in (c) 
indicate the edge states are chiral. (e) A weaker hopping be- 
tween the edge (dashed black vertical bonds) is reintroduced, 
(f) The gapped energy spectrum corresponds to (e). (g) An 
edge soliton (marked by the gray plaquette) is introduced by 
reversing the sign for half of the vertical bonds (red dashed 
lines) between the edges, (h) The energy spectrum corre- 
sponds to (f). The spectra in (b), (d) and (h) are obtained 
from H F with Ji = J 2 = 1, J3 = 0.2 and J 4 = 0.5. The 
spectrum in (f) and (h) was obtained with a restored edge 
coupling Jz — ±0.02. 



by a massive Dirac Hamiltonian 

H E = J dx (~iv^ a z d x ip + mi^ a x ip), (8) 

where v is the edge velocity, xp^ = (V']{!'0l) with 
being the right /left fermion creation operators. When 
the restored bonds have a sign reversal (Fig. [3g), the 
mass term in Eq. ([8]) becomes x-dependent and changes 
sign as x goes through the location of the topological 
defect. This should result in one localized zero mode per 
defect according to Ref.@, Fig. ([3ji) shows that this 
is indeed true. The presence of such zero mode plus the 
fact that the spectra of Majorana fermion models are 
E «-> —E symmetric, immediately imply that the zero 
modes are Majorana fermion states in the original Kitaev 
model. Furthermore, using the argument of Ref. [24j it 
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FIG. 4: The evolution of the of the ft-surface as the model 
in Fig. (Tib) is gradually deformed into Haldane's model. In 
constructing the figure we used Ji = J2 = J3 = 1. The full 
strength of the second neighbor hopping is 0.5. 



can be shown that the braiding of such Majorana fermion 
defects leads to the non-Abelian statistics. Thus, via 
the mechanism of Refs.jH, HJ a two dimensional defect 
with fractionalized quantum number has emerged! Its 
presence is determined by the topological nature of the 
host bulk band insulator just as the edge states are. 

Actually, a similar phenomenon was also found in the 
continuum theory of Dirac fermions interacting with the 
topological defects of a Higgs field[25[. However, in all 
the examples we considered here, the location of the edge 
Dirac point in the momentum space is far away from 
those of the bulk Dirac points. Consequently the theory 
discussed in Ref. is not applicable here. 

As a digression, we now show the fermion represen- 
tative of the Kitaev model is topologically equivalent 
to Haldane's lattice model for integer quantum Hall 
effect [l2j . Fig. Q shows the evolution of the /i-surface 
by gradually switching off the second neighbor hopping 
in Haldane's model which are not contained in Eq. ((6|. 
The leftmost column are the /i-surfaces for the Haldane 
model (top) and Kitaev (bottom), respectively. In the 
rest of the figure the surfaces are dissected to reveal the 
origin (the red dot). Clearly as we follow the evolution 
the origin never migrate across the /i-surface. Thus V for 
the two models are the same. 

Finally what about the Laughlin quasiparticles? Al- 
though fractional quantum Hall liquids are not band in- 
sulators they are clearly topological insulators. Indeed, 
as shown by Wen [2r3 |. when a quantum Hall liquid on a 
torus is cut open (Fig. [3k), there are "chiral Luttinger liq- 
uid" in-gap edge modes. At 1/m filling, the edge modes 
are described by the following free boson Hamiltonian 



h b = J dx{|ri(£) 2 + £[cw>(*)] 2 }, 



tial V — — g J dx cos(m(f>) needs to be added j26|. When 
g is sufficiently big, a gap opens in the edge spectrum 
and the ground states become m-fold degenerate. They 
are characterized by (4>(x)} = 2-Kl/m{l — 0, ...,m — 1). 
In this case an edge soliton is where <fi(x) interpolates 
between two different ground states, say, (<j){x)) = and 
2ir/m. As shown by Goldstone and Wilczek 27|, such a 
soliton carries a charge AQ = A<fi/2ir = l/m,precisely 
the same as a that of a Laughlin quasiparticle pj . Hence 
the Laughlin quasiparticles are also edge solitons of a 
topological insulator! 
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(9) 



where n and <j) are conjugate boson fields satisfying 
[n(x), 4>(y)) — i5(x — y). To reconnect the edges, a poten- 



[1] 
[2] 
[3] 

[4 
[5 
[6 

[T 
[8 
[9 
[10 

[11 
[12 

[13 

[14 

[15 
[16 

[17; 

[is; 

[19 

[20 

[21 

[22 
[23 
[24 
[25 

[26 
[27 



R. B. Laughlin, Phys. Rev. Lett. 50, 1395 (1983). 
R. Jackiw and C. Rebbi, Phys. Rev. D 13, 3398 (1976). 
W. P. Su, J. R. Schrieffer and A. J. Heeger, Phys. Rev. 
B 22, 2099 (1980). 

F. D. M. Haldane, Phys. Rev. Lett. 66,1529 (1991). 

G. Moore and N. Read, Nucl. Phys. B 360, 362 (1991). 
N. B. Kopnin and M. M. Salomaa, Phys. Rev. B 44, 9667 
(1991). 

G. E. Volovik, JETP Lett. 70, 609 (1999). 
N. Read and D. Green, Phys. Rev. B 61, 10267 (2000). 
M. Stone and R. Roy, Phys. Rev. B 69, 184511 (2004). 
S. Tewari, S. Das Sarma and D.-H. Lee, Phys. Rev. Lett. 
99, 037001 (2007). 

A. Kitaev, Ann. Phys. 321, 2 (2006). 

F. D. M. Haldane, Phys. Rev. Lett. 61, 2015 (1988). 

C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, 226801 

(2005) . 

S. Murakami, N. Nagaosa and S.C. Zhang, Science 301, 
1348 (2003). 

A. Seidel et al, Phys. Rev. Lett. 95, 266405 (2005). 

A. Seidel and D.-H. Lee, Phys. Rev. Lett. 97, 056804 

(2006) . 

E. J. Bergholtz and A. Karlhede, Phys. Rev. Lett. 94, 
026802 (2005). 

E. J. Bergholtz et al, Phys. Rev. B 74, 81308 (2006). 
X.-Y. Feng, G.-M. Zhang and T. Xiang, Phys. Rev. Lett. 
98, 087204 (2007). 

G. Baskaran, S. Mandal, and R. Shankar, Phys. Rev. 
Lett. 98, 247201 (2007). 

W.-Y. Hsiang and D.-H. Lee, Phys. Rev. A 64, 052101 
(2001). 

D. J. Thouless et al, Phys. Rev. Lett. 49, 405-408 (1982). 
Y. Hatsugai, Phys. Rev. Lett. 71, 3697-3700 (1993). 

D. A. Ivanov, Phys. Rev. Lett. 86, 268 (2001). 

C. G. Callan Jr. and J. A. Harvey, Nucl. Phys. B 250, 

427 (1985). 

X. G. Wen, Phys. Rev. B41, 12838-12844 (1990). 

J. Goldstone and F. Wilczek, Phys. Rev. Lett. 47, 986- 

989 (1988). 



